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We study the magnetocaloric effect for the integrable antiferromagnetic high-spin chain. We 
present an exact computation of the Griineisen parameter, which is closely related to the magne- 
tocaloric effect, for the quantum spin-s chain on the thermodynamical limit by means of Bethe 
ansatz techniques and the quantum transfer matrix approach. We have also calculated the entropy 
S and the isentropes in the (H, T) plane. We have been able to identify the quantum critical 
points He = — rr looking at the isentropes and/or the characteristic behaviour of the Griineisen 

s + 2 

parameter. 
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I. INTRODUCTION 

The magnetocaloric effect has been known for more 
than hundred yearsi and it is related to the temperature 
change of magnetic systems induced by an adiabatic vari- 
ation of the external magnetic field. In recent years it has 
received considerable attention in view of new potential 
cooling applications^. 

Many different families of magnetic materials ranging 
from ferromagnetic, ferrimagnetic to antiferromagnetic 
systems have been shown to present large or unusual 
magnetocaloric effect 3 ". In particular, the magnetocaloric 
effect has been measured in (quasi) one-dimensional ma- 
terials which behave as quantum spin- 1/2 chains 4 - and 
high-spin chains^. 

Moreover the existence of the magnetocaloric effect in 
onc-dimcnsional systems has been studied theoretically 
for spin-1/2 Heisenberg chain£, for spin-s chaini and for 
mixed-spin chains^ by means of numerical calculations. 
There are also some exact results for the XY chain in 
transverse field^ and for the Ising model^. Recently the 
magnetocaloric effect and the isentropes in the magnetic 
field/temperature (H,T) plane have been obtained ex- 
actly by Bethe ansatz techniques for the integrable spin- 
1/2 Heisenberg chain 1 ^. 



The magnetocaloric effect (§jj) s = — ^J/dT) an< ^ 
the related quantity called Griineisen parameter T h have 
been pointed out as an important tool to detect and clas- 
sify quantum critical points&ii. The Griineisen parame- 
ter for a magnetic systems can be written 

Th = ^{m) s = -^{w) h ' (1) 

where Ch is the specific heat at a constant magnetic field 
and (^f-)jj is the temperature variation of the magneti- 
zation M. This parameter Th has a characteristic sign 
change close to the quantum critical point, which is due 
to the accumulation of entropy at the critical points. 

The integrable spin-s generalization of the Heisen- 
berg models was exactly solved long agoi 3 - providing 
all the eigenvalues and eigenvectors in terms of the 
Bethe equations. Its thermodynamic properties have 
been firstly studie d 13 ' 14 by means of the thermodynamic 
Bethe ansatz (TBA) method-^, which consist of an infi- 
nite number of non-linear integral equations (NLIE) for 
the free-energy. Alternatively, using the quantum trans- 
fer matrix (QTM) approach^, it was derived a finite 
number of NLIEil which is more suitable for practical 
calculations. 

Here we are interested in the exact computation of the 
magnetocaloric effect (Th), entropy and the isentropes 
in the (H, T) plane for the integrable antiferromagnetic 
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spin-s chain. The thermodynamic quantities required to 
achieve this goal, like entropy, specific heat and magneti- 
zation, are determined as a function of temperature and 
magnetic field by means of the solution of a finite set of 
NLIE which arises from the QTM approach^. 

This paper is organized as follows. In section [TTl we 
outline the integrablc Hamiltonians and the associated 
integral equations. In section ITO1 we present our results 
for the magnetocaloric effect and the isentropes in the 
(H, T) plane. Our conclusions are given in section IIVI 

II. HAMILTONIAN AND INTEGRAL 
EQUATIONS 

The Hamiltonian of the integrable spin-s generaliza- 
tion of the Hcisenbcrg model for s = 1/2, 1 and 3/2 are 
given by 
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where L is the number of sites and Si = (Sf , Sf , Sf) are 
the SU{2) generators. 

One can write a closed formula for the Hamiltonian 
assuming a generic spin-s value a follow, 



(5) 



where 



Qs(x) = J2 [2V(J + 1) - V-(l) - V>(2s + 1)] I 



3=0 



k = 

k^3 



X Xfc 
Xj Xfc 



(6) 

with Xk — \ [k(k + 1) — 2s(s + 1)], ip(x) the digamma 
function and J is the exchange constant. From now on we 
assume J = 1. Note that we have also added a Zeeman 
term on the Hamiltonian ([5]). 

The free-energy of the system per lattice site calculated 
at the thermodynamic limit (L — > oo) is given by 



f(T,H) = f -T(K*hiBB) (0), 



(7) 



where f = i>(2s + 1) - V(^) + V>(|) - ^(1), if (a?) 
coSK\nx\ anc ^ ^ e s y m bol * denotes convolution / * g(x) 

(4) JZo f( x ~ v)a(y)dy. 



The auxiliary functions b(x), b[x) and its simply related functions B(x) = b(x) + 1 and B(x) = b(x) + 1 are solution 
of the following set of non-linear integral equations^ 



lny( s ~^(x) 
]nb(x) 
In b(x 



-)8d(x) + /3f 
y -/3d(x-) - 0§ J 



( bxY&(x) \ 

lnB(x) 
V lnB(x) 



(8) 



where = 2 cosh [nx] » ^ = i s the inverse of temperature and i? is the magnetic field. 
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The kernel matrix is given explicitly by 



K{x) 



\ 





K{x) 
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K{x) 





which is a matrix of dimension (2s + 1) x (2s + 1) with F(x) 





K{x) 


K(x) 
K{x) 

= 11 






K(x) 
F(x) 
-F(x - i) 

-\k\/2+ih 






K(x) 
-F(x + i) 
F(x) J 



(9) 



dk. 



-oo 2 cosh [k/ 2] 

In order to obtain the desired thermodynamical quantities, we can calculate the derivatives of the free-energy with 
respect to temperature T (or more conveniently /3) and magnetic field H. It turns out to be more efficient to calculate 
the derivatives of the free-energy in terms of the solution of linear integral equations. These equation are obtained by 
differentiation of the equation (|8]). This way we can avoid numerical differentiation of the free-energy. 

Specifically one can write the entropy S — — (^§j^j = P 2 (f^) as f°ll° w J 

S = {K *lnBB)(0) - /3(K * dp\nBB)(Q), (10) 



where dphxB(x) = ^^[dp hib(x)}, dplnB(x) - 
j = 1/2, . . . , s — 1/2 . These new auxiliary functions are solution of the following system of linear integral equations 

/ dplny^\x) \ ( \ / dp\nY^{x) \ 



b(x) 



T [Spin 6(a)] and dp]nY^(a 



^[dp]ny^(x)] for 



dp lny( s ~^(x) 
dp ln6(x) 
^ dp\nb(x) J 



-d(x) + f 
y -d(x) - § J 

To obtain the entropy in the (H, T) plane, one has to solve the above equations ((8]) and (fTTj) varying the temperature 
and the magnetic field. 



dpkiY^-i^x) 
dphiB(x) 
\ dp\nB{x) J 



(11) 



The specific heat C H =T (|f ) H 



-j8 



can be obtained from (fTOI) . 



C H = /3 2 (K*d 2 p\nBB){0), 
which is given in terms of the solution of following linear integral equations 



(12) 



/ a§lny(i)(x) \ 



dglny( s s)(x) 
dpnb(x) 
\ dpnb(x) ) 



= K* 



( dj]nY(i\x) \ 

d 2 p \nY^-i\x) 
dj\nB{x) 
\ djlnB(x) J 



(13) 



where 



^B(x) = ^ {^g^ + [d}lnb(x)]} and flg hrWfr) = { + [agln^)( g )]}. 



In order to obtain the Griineisen parameter we have also to determine the {j^r) H - Therefore we have firstly 
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to calculate the magnetization M = 

M = j(K*d H \nBB)(0), 



|£) from©, 



(14) 



which is written in terms of the auxiliary function 
8h^B(x) = [9h lnb(x)] (likewise for the other 

auxiliary functions) that are now derivatives with respect 
with the magnetic field. These new auxiliary functions 
are solution of the following system of linear integral 
equations 

\ / n \ / 



/ 



8h In y' 2 ' (x) 



8h lny( s ~^)(x) 
&h In b(x) 
y d H In b(x) J 



-K* 



d H ]nB(x) 



\ d H lnB(x) J 
(15) 

The derivative of the magnetization with respect 
to temperature for constant magnetic field {^jr) H = 
—(3 2 (j^jfj can be finally obtained from (jT4j) , which 
results 

!^L=(K* d H lnBB)(0) - 0(K * d} H ]nBB)(0), (16) 



where &i H In B(x 



b(x 



PH iix " VV — b(x) 

which should satisfies 
( dj H ]ny&(x) \ ( n \ 



lnb{x)][d H lnb(x)] 
b(x)+l 



[dj H \nb{x)]) 



(■v) 

^ dj H \nb{x) J 



dj H ]nb(x) 

32 



-K* 



( d} H \nY^\x) \ 

dj H lnY^ s -^(x) 
8} H \nB(x) 

(17) 



III. GRUNEISEN PARAMETER AND 
ENTROPY 
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In this section we will present the results for the 
Griineisen parameter, which is closely related to the mag- 
netocaloric effect. We will also show the results for the 
entropy and the isentropes in the (H, T) plane. 

First we show on Fig. ([T]) the Griineisen parameter 
Th as a function of magnetic field for some spin values 



FIG. 1: (Color online) The Griineisen parameter Fh for dif- 
ferent values of the spin for (a) T = 0.05, (b) T — 0.5 and (c) 
T = 5. 

s = 1/2, 1,3/2, 2 and 5/2. Note that the case s = 1/2 was 
calculated on Ref. 10. For the temperature T — 0.05, we 
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FIG. 2: (Color online) Entropy S(H,T) for s = 1,3/2. The 
isentropes are for 5* = 0.1, 0.2, . . . , 1.3. 

can sec that the transition to saturation at He = —*rr 

S+ 2 

is signaled by sign changes of the Griineisen parameter 
from negative to positive values toward the higher fields 
values. For higher temperatures, like T = 0.5, we note 
that these sign changes move away from the zero temper- 
ature saturation field, which separates the antifcrromag- 
netic and ferromagnetic phases. Finally, if we go further 
to higher temperatures, e.g T = 5 one can see that all the 
characteristic behaviour have disappeared, which imply 
that the thermal fluctuations are already strong enough 
to drive the system to excited states where no quantum 
phase transition effects can be seen. Moreover there is a 
small structure at low magnetic fields and low tempera- 
tures Fig. HJa) which is due to the singular nature of the 
point (H = 0) of the isotropic integrable spin chains^. 



FIG. 3: (Color online) Entropy S(H, T) for s = 2, 5/2. The 
isentropes are for 5* = 0.1, 0.2, . . . , 1.6. 

We show on Figs. [5] and [3] the entropy and the 
isentropes for the spin-s chain in the (H, T) plane for 
s = 1,3/2,2 and 5/2. 

The quantum phase transitions are indicated by the 
isentropes which are tilted towards the quantum critical 
points showing a minima nearby He or equivalently the 
entropy peaks at the critical point. This accumulation 
of entropy nearby the critical point indicates that the 
systems is maximally undecided which ground state to 
choose^. Moreover the Griineisen parameter, which is 
proportional to the slope of the isentropes (f§0, has a 
different sign on each side of the quantum critical point 
as we have shown on Fig. [TJ Besides that the isentropes 
are very steep nearby the critical point indicating the 
existence of a large magnetocaloric effect. 
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IV. CONCLUSION 

In this paper we have studied the magnetocaloric effect 
for the integrable spin-s chain. We have calculated the 
Griincisen parameter, which is proportional to the mag- 
netocaloric effect, as a function of the external magnetic 
field on the thermodynamic limit and at finite tempera- 
tures. We have also obtained entropy and the isentropes 
in the (H, T) plane. 

( s) 

The quantum critical point He have been identified 
by the minima of the isentropes and by the sign changes 
of the Gruneisen parameter as a function of the magnetic 



field. Our results are in agreement with the previous 
results for the s = 1/2 case^. 

We hope that our exact results could be useful 
for understanding experimental results for quasi one- 
dimensional systems, e.g^. We also expect that our re- 
sults could be further extended to the case of alternating 
spin- (Si, 52 ) chami£. 
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